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A status report of utilizing soft electromagnetic radiation (aka ther- 
mal photons and dileptons) in the diagnosis of strongly interacting matter 
in ultrarelativistic heavy-ion collisions is given. After briefly elaborating 
on relations of the electromagnetic spectral function to chiral symmetry 
restoration and the transition from hadron to quark degrees of freedom, 
various calculations of electromagnetic emission rates in the hadronic and 
quark-gluon plasma phases of QCD matter are discussed. This, in particu- 
lar, includes insights from recent thermal lattice QCD computations. Ap- 
plications to dilepton and photon spectra in heavy-ion collisions highlight 
their role as a spectro-, thermo-, baro- and chrono- meter of extraordinary 
precision. 

PACS numbers: 12.38. Mh,21.65.Jk,25. 75. Cj,25.75.Nq 

1. Introduction 

Collisions of heavy nuclei at (ultra-) relativistic energies provide the fas- 
cinating opportunity to recreate blobs of strongly interacting matter which 
last existed naturally almost 14 billion years ago, for a few microseconds af- 
ter the Big Bang. They are furthermore the only means by which the bulk 
properties of a nonabelian gauge theory can be compared to experiment. 
This is particular important in a situation where the coupling strength is 
not small, implying a wealth of nonperturbative and collective phenomena 
which can only be unraveled in close collaboration between experiment and 
theory. 

Over the last twenty years or so, ultrarelativistic heavy-ion collisions 
(URHICs) at various laboratories around the world have demonstrated that 
systematic studies over a large region of the QCD phase diagram are possi- 
ble, cf. left panel of Fig. [TJ For example, the measured hadron abundances 



Presented at 51. Cracow School of Theoretical Physics on "Soft Side of the LHC" 



(1) 



2 



ZAKOll PRINTED ON OCTOBER 20, 2011 



0.25 
0.20 
> 0.15 

CD 

C5 

H 0.10 
0.05 
0.00 



Heavy Ions 



QGP 

quasiparticles, 
resonances? 

<qq>=<qq>=0 



HG ^ 

7i,p,NA... r CSC 

<qq>+0 \Coopei>ijairs,.. 

1 1 <qq>*o 



u, d, s 



- 3-loop pQCD 
Naive quark model 



y5(GoV) 



0.0 0.5 1.0 1.5 2.0 

[GeV] 

Fig. 1. Left panel: schematic phase diagram of strong- interaction matter, indicating 
a cross-over transition (dashed hne) at large temperature, T, and smaU baryon 
chemical potential, fiB, and a first order transition into a color-superconductor 
(CSC) at small T and large iib (solid line). The "data points" represent {fiB,T) 
values extracted from ratios of different hadron species employing a thermal hadron- 
resonance-gas model [T]. Right panel: EM spectral function in vacuum as measured 
by the i?-ratio in e~^e~ annihilation into hadrons (figure taken from Ref. [2])- The 
EM spectral function is the only one whose medium modifications are directly 
accessible in heavy-ion collisions, through dilepton invariant-mass spectra. 



can be well explained by a chemically equilibrated hadron-resonance gas in- 
cluding strange particles (contrary to p-p collisions) [IJ, and the transverse- 
momentum (pt) spectra of hadrons can be well described by an explosive 
blast-wave source with common temperature and collective "flow" velocity 
in excess of half the speed of light. The observed azimuthal asymmetries in 
the Pt spectra in non-central collisions (the so-called elliptic flow, V2) can 
be well accounted for through the pressure gradients as obtained from a lo- 
cally equilibrated medium driven by hydrodynamic expansion [3]. However, 
hadronic observables ultimately emanate from the "freezeout" configuration 
of the fireball (typically at a temperature of Tfo ~ 100 MeV), and thus do 
not provide information on the microscopic structure and interactions of 
medium. 

Electromagnetic (EM) observables, photons and dileptons, are special 
(see Refs. [4,^ for recent reviews). Their interaction rate in the strongly 
interacting medium is small enough for them to escape the interior of the 
fireball unaffected, but large enough to be produced in measurable quan- 
tities. The radiation of photons from a thermalized fireball has been long 
recognized as powerful thermometer. The emission of virtual photons (aka 
dileptons, e~^e~ or ii'^p,^) carries additional information in terms of their 
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invariant mass, M = {q^ — (fY^"^. In fact, dilepton invariant-mass spectra 
are the only observable which gives direct access to the in-medium mod- 
ification of a hadronic spectral function (strong decays like A — t- ttN or 
p — )• TTTT suffer from final state absorption and are largely emanating from 
the freezeout process of the collision) . This is evident from the 8-differential 
thermal production rate, 

"^^'^ - "™ /^(<Zo;T)ImnEM(M,g;/XB,r), (1) 



which only depends on the in-medium spectral function, Im IIemj its thermal 
weight in terms of the Bose-Einstein distribution, and a free virtual- 
photon propagator, A clean temperature measurement becomes pos- 
sible if ImllEM is reliably known, which is usually the case for large 
where temperature corrections can be assessed perturbatively (in principle 
this also applies for photons at large = g^, but the leading term is al- 
ready nontrivial in the strong coupling, i.e., ©(a^T^), while for dileptons it 
is 0(1)). However, the physics accessible through ImllEM is much richer. In 
the vacuum, the latter is accurately known from the inverse process of e~^e~ 
annihilation into hadrons, as displayed in the right panel of Fig.[T]in terms of 
the famous ratio R = (T(e^e~— 7>hadrons)/(T(e"'"e~— )-/x"''//^) oc Imllg^/M^. 
Its medium modifications encode a wide variety of properties of the strongly 
interacting medium. For example, transport coefficients can be extracted 
by the spacelike and timelike limits of vanishing energy and momentum, 
corresponding to the EM susceptibility and conductivity, respectively, of 
the medium. A more ambitious goal is to systematically map out how the 
hadronic degrees of freedom, represented by the p, uj and resonance peaks, 
are affected across the QCD phase diagram and in this way reflect its struc- 
ture. In particular, one would like to find out how the nonperturbatively 
generated mass of the vector mesons dissolves and how they ultimately yield 
to a spectral function of weakly correlated quarks at high temperatures, sig- 
naling deconfinement. Note that light vector mesons (especially the p) may 
well serve as a prototype hadronic resonance (it is neither a Goldstone boson 
nor a heavy-quark state) in the sense that their fate in hot/dense matter is 
shared by large set of resonances. A hadron-resonance gas is known to give 
a good approximation of the the QCD partition function until close to the 
phase transition [6]. 

Mass de-generation and deconfinement are closely related to the fate 
of the QCD condensates in the medium. Figure [2] shows recent results 
of thermal lattice QCD for order parameters of the corresponding QCD 
phase transition(s). These calculations are carried out for A'^y=2+l-flavor 
QCD with realistic light quark masses and clearly demonstrate the cross- 
over of the transition from a chirally broken hadronic phase to chirally 
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Fig. 2. Results of lattice-QCD computations for the renormalized Polyakov 
loop (left) and for the "strangeness-subtracted" chiral condensate normalized to 
its vacuum value, {qq){T)/ {qq){0) (right). Both quantities are order parameters 
of the QCD phase transition, the former associated with quark deconfinement in 
the limit of large quark masses and the latter with chiral symmetry restoration 
for ruq — > 0. The computations are carried out with realistic quark masses for 
Nf = 2 + 1 flavors with the grey band indicating the extrapolated continuum limit. 
Figures are taken from Ref. [TJ [S] . 

restored partonic matter. It is rather intriguing, though, that the inflec- 
tion points of these quantities, usually associated with the pseudo-critical 
temperature, seem to be separated by about 20 MeV, i.e., ~ 150 MeV 
and T™°^ ~ 170 MeV for chiral restoration and deconfinement, respectively. 
Even without the notion of a pseudo-critical temperature, which is not un- 
ambiguous, it is still remarkable that at, say, r~160 MeV the Polyakov 
loop has risen by merely ~10% toward its asymptotic value while the chi- 
ral condensate has already dropped to below 50%. This is suggestive for a 
medium in which chiral symmetry is largely restored but hadrons are still 
prevalent degrees of freedom (not unlike the "quarkyonic phase" conjectured 
to exist in the moderate- /i^ region of the QCD phase diagram [9j). This 
would furthermore imply that hadronic degrees of freedom are the adequate 
basis to address mechanisms of chiral restoration. Partial chiral restoration 
in the hadronic phase is very encouraging from an experimental point of 
view, as it turns out that thermal dilepton emission in the low-mass region 
(LMR, M < 1 GeV) is dominated by the contribution from hot and dense 
matter at temperatures T ~ 150 — 200 MeV (this follows from the interplay 
of the increasing 3-volume and the decreasing thermal Bose factor in the 
space-time integrated spectral yield, see Sec. l3.1l below) [10]. 

The manifestation of the Spontaneous Breaking of Chiral Symmetry 
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Fig. 3. Vacuum spectral functions of the vector and axialvector currents as mea- 
sured in hadronic r decays into Goldstone bosons [131 [T3]. The hues are model 
calculations using p and ai dominance with smoothly onsetting multi-pion con- 
tinua. 

(SBCS) in the vacuum EM spectral function becomes apparent upon de- 
composing it into good isospin states. The dominant contribution arises 
from the isovector channel (/=1) which in the LMR is essentially saturated 
by the p meson (the contribution of the isoscalar-vector (a;) channel to Hem 
is down by a factor of ~10). The Weinberg sum rules |11[ [T2] relate mo- 
ments of the isovector-vector and -axialvector spectral functions to order 
parameters of SBCS, 

oo 

fn = - f-s^ [ImUvis) - ImUAis)] . (2) 
J vr 



Both spectral functions have been accurately measured in vacuum via hadro- 
nic decays of the r lepton at LEP [131 E] (see Fig. [3|) and constitute one 
of the best experimental evidences for SBCS. For example, for n = — 1, 
one has /_i = where /7r=93MeV is the pion decay constant (or pion 
"polestrength" ) . Recalling the Gellmann-Oakes-Renner relation, /^m.^ = 
— mg(0|^g|0), one recognizes the close relation between scalar quark conden- 
sate, Goldstone bosons and the isovector axial-/vector spectral functions 
(the explicit chiral breaking is signified by the small current quark mass, 
niq ~ 5MeV, and the pion mass). The Weinberg sum rules, which remain 
valid in the medium, exhibit the important role of the axialvector spectral 
function in the search for chiral restoration from dilepton data: the vector 
spectral function degenerates with its chiral partner, the axialvector, in the 



6 



ZAKOll PRINTED ON OCTOBER 20, 2011 



chirally restored phase. Similar to the vector channel, the nonperturbative 
part of the axialvector channel is dominated by a resonance, the ai(1260). 
At higher masses, M > 1.5 GeV, one expects from perturbation theory that 
the axialvector also merges into a continuum which coincides with the vector 
channel. Another example of relating QCD condensates to empirical spec- 
tral are QCD sum rules whose application to the p meson will be discussed 
below. 

The remainder of this article is organized as follows. In Sec. [2] we intro- 
duce the vector-current correlator in strongly interacting matter and recall 
its non-interacting and vacuum limits (Sec. 12. ip . review calculations of its 
spectral function in hadronic (Sec. 12. 2p and partonic media including recent 
results from thermal lattice QCD (Sec. 12.3]) . and discuss the corresponding 
thermal production rates of dileptons and photons (Sec. 12.4]) . In Sec. [3] we 
turn to the quantitative analysis of EM emission spectra in URHICs, be- 
ginning with descriptions of the space-time evolution of the medium over 
which the rates need to be integrated (Sec. 13. ip . Theoretical calculations 
of EM spectra in URHICs will be compared to experiment to extract in- 
formation on medium effects of the spectral function, the time duration of 
emission (spectro- and chrono-meter, Sec. 13. 2p and temperature and collec- 
tive properties of the expanding fireball (thermo- and baro-meter, Sec. 13. 3p . 
We conclude in Sec. [H 



2. Vector-Current Correlator in Medium 

In this section the central role is played by the hadronic correlation 
function of two electromagnetic currents defined as 

^EuiQ) = -^l d'x e^'^^ e(xo) ([j^M(^)>iEM(0)]>T , (3) 

where (• • ■)t denotes the expectation value at finite temperature (giving 
rise to a retarded spectral function). In a partonic basis, corresponding to 
the elementary degrees of freedom in the QCD Lagrangian, the EM current 
takes the form 

Jem = ^^^""u - -d^d - -s-i^s (4) 

where we constrain ourselves to the 3 light flavors up, down and strange. 
The EM current in quark basis, Eq. (jj]), can be rearranged into good isospin 
states which naturally leads to the hadronic basis according to 

Jem = \{u-i^'u-d-f^'d) + ^{u-i^'u + d-i''d)-^sYs 
= -I- ^ 1^ - - i^ 
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with the properly normahzed hadronic currents, {v = p,U},(f)). Invok- 
ing the vector-dominance model (VDM) characterized by the field-current 
identity = m^/gy, the hadronic currents are saturated by the light vec- 
tor mesons and the pertinent current-current correlator turns into a vector- 
meson propagator, D^. 



2.1. EM Spectral Function in Vacuum and Thermal Emission Rates 

The imaginary part of the EM current-current correlator is the EM 
spectral function. In the vacuum it is directly observable in e'^e~ anni- 
hilation, recall right panel of Fig. [TJ Its structure suggests a decomposi- 
tion into a nonperturbative low-mass regime, M < 1 GeV, saturated by 
hadronic resonances, a "dip" region for 1 GeV< M < 1.5 GeV, and regime 
for M > 1.5 GeV where perturbation theory applies, 

_E (^)' ImZ)r(M) ,M<Mrai, 



ImnSS(M) 



=p,UJ,(f> 



q=u,d,s 



(6) 

{Nc=3 is the number of colors and Cq the quark charge in units of the 
electron charge). Note that in the perturbative regime the final state in e^e~ 
annihilation is still hadronic (with broad overlapping excited resonances 
such as p', uj', etc.), but the "long-distance" hadronization process does 
not affect the short-distance quark-production cross section. This is often 
referred to as parton-hadron duality, which in the vector channel is realized 
for masses beyond a "duality threshold", M > MJjJ^j ~ 1.5 GeV. 

The inverse process of e~e^ annihilation is intimately related to the 
production rate of thermal dileptons from a heat bath of strongly interacting 
matter. It can be calculated starting from Fermi's Golden rule. 

Fee dNee 4 /" TT ^^P^ TT ^^Pf n Xf D JD \ 

x(^iiri/)(/brN)r(i ± /^) (ie^^j:) > (7) 

which can be cast into the form given in Eq. ([1]). This expression is leading 
order in the fine structure constant, 0{a'^^), but exact in the strong inter- 
actions which are encoded in the spectral function (the leading order in the 
strong interaction is 0(1)). In the same framework one can calculate the 
thermal production rate of real photons as 

10 = f''{qo;T) ImUly^iM = 0,q = qo;PB,T) , (8) 
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which only involves the transverse polarization of the EM spectral function. 
The leading orders of the photon rate are one lower in the EM coupling, 
0{aEM), but one higher in the strong coupling, 0{as), relative to the dilep- 
ton rate. This implies significant uncertainties in the magnitude of the rate 
even in the perturbative regime. 

In the following three sections we focus on the dilepton rate in the soft 
(low-mass) regime, where the key issue is to determine the in-medium EM 
spectral function. 

2.2. Hadronic Matter 

At low temperatures and densities, the appropriate degrees of freedom 
to describe strongly interacting matter are hadrons. In lukewarm matter, 
T < m^, at vanishing baryon chemical potential, ^s, the most abundant 
species are pions whose low-energy interactions are determined by the chiral 
Lagrangian. At higher temperatures, T > mj^, hadronic resonances quickly 
outnumber pions; much less is known about the chiral structure of resonance 
interactions (including excited nucleon states). However, it is still possible 
to write effective hadronic Lagrangians compatible with basic symmetry re- 
quirements and constrain their parameters by empirical information, e.g., 
hadronic and radiative decay widths. The such constrained interactions 
can then be implemented into many-body theory to calculate in-medium 
selfenergies and spectral functions of hadronic excitations. Clearly, the re- 
liability of the predicted spectral functions largely hinges on the quality of 
the constraints applied to the effective interactions. 

The above strategy has been adopted by several groups to compute 
the in-medium spectral function of the p meson (cf., e.g., the review pa- 
pers [151 HSl HI El US] for a detailed discussion and further references). 
The special role of the p meson simply derives from its dominant con- 
tribution to the low-mass EM correlators where it outshines the cj by a 
factor of ~10 (the constituent-quark model predicts a factor of 9, as in- 
ferred from Eq. ([5]); from the empirical dilepton decay widths one finds 
= ("T'p/5p)/("iw/fi'S) —11)- Starting point is a realistic de- 
scription of the p properties in the vacuum. Usually the p is introduced 
into the chiral Lagrangian via a local gauge symmetry. To one loop order, 
the vacuum p selfenergy is built up from the standard p — )• vrvr diagram and 
a tadpole arising from the ppmr 4-point vertex (the pertinent diagram has 
no imaginary part). With 3 free parameters (bare p mass, gauge coupling 
g and a regulator (cutoff) for the divergent loop integrals), a satisfactory 
reproduction of the p spectral shape as seen in the pion-EM formfactor and 
P-wave TTTT scattering can be achieved, cf. Fig. U] (the r — )• 2^1^^ data are also 
well described, see Fig. A good description of the entire spectral shape 
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Fig. 4. Electromagnetic formfactor of the pion (left panel) and P-wave tttt scatter- 
ing phase shifts (right panel) in the vacuum. An effective model assuming p-meson 
dominance with coupling to two-pion states is compared to data. Figures taken 
from Ref. [19]. 



- not just its mass and width - is essential for later purposes of calculating 
low-mass dilepton spectra, especially below the free p mass. 
Medium modifications of the p propagator, 

Dp{M, q- pB,T) = —2 ^— ^ — , (9) 

are induced by (a) interactions of its pion cloud with hadrons from the heat 
bath, included in SpTTTr (e.g., 'rN — A), and (b) direct p scattering off 
mesons (e.g., pvr — )• ai) and baryons (e.g., pN — )• A^*), denoted by Y^pu and 
YjpB , respectively. The medium effects in nuclear matter can be comprehen- 
sively constrained by analyzing vrA^ — >• pN scattering and photoabsorption 
spectra on the nucleon and nuclei [20]. In this context, the pion cloud effects 
correspond to meson-exchange currents (e.g., pion exchange in — )• 7A), 
and the direct nucleon excitations, 7A'^ — )• A'^*,A*, produce the resonance 
peaks in the cross section. For the p interactions in a meson gas, one has 
to resort to the hadronic and radiative decay branchings of mesonic res- 
onance with large coupling to p and 7 final states, e.g., oi — )• vrp, 7r7 or 
K\ — >■ Kp^K^ [2IJ. To date, for the case of cold nuclear matter, different 
calculations of the in-medium p spectral functions have reached agreement 
at a semi-quantitative level, and remaining discrepancies can be largely 
traced back to slight variations in the treatment of medium effects (e.g., 
missing pion-cloud effects or linear- in-density approximations). An example 
for the p spectral functions under hot and dense conditions as expected for 
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Fig. 5. Spectral function of the p meson in hot hadronic matter; left panel: in hot 
and dense matter for temperatures r=120, 150 and 180 MeV and at fixed baryon 
chemical potential of /^s =330 MeV, corresponding to baryon densities of O.lgoj 
Q.Tqq and respectively (£io=0.16fm~'^); right panel: in a hot meson with all 

baryon-induced medium effects switched off. 

heavy-ion collisions at the CERN-SPS is displayed in the left panel of Fig. [5j 
The p resonance peak undergoes a strong broadening, indicative for its ul- 
timate "melting" close to the phase transition. The medium modifications, 
especially the low-mass enhancement, are much reduced if baryon-induced 
effects are switched off, see right panel of Fig. [5l 

Let us briefly allude to medium effects in the dip region, i.e., for masses 
between 1 and 1.5 GeV. In this regime the continuum starts to develop, 
characterized by multi-meson contributions in the spectral function whose 
interactions with a medium are difficult to assess microscopically. Fortu- 
nately, one can make a more simple yet elegant argument based on chiral 
symmetry to estimate the medium effects on the spectral function. It was 
first developed in Ref. [22] for the finite-temperature case. Using current 
algebra in the chiral limit (771,7^ = 0), the interactions of the vector and axi- 
alvector correlators with a lukewarm pion gas were shown to result in their 
mutual mixing as 

Ilv{q) = {l-e) liTiq) + e IiT{q) (10) 

(and likewise for the axialvector upon exchanging V ■v^ A), with the mixing 
parameter e=r^/6/^. The model-independent leading-temperature effect 
on the vector spectral function is an interaction with a thermal pion which 
reduces the strength at its resonance and moves it into the axialvector res- 
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Fig. 6. Effect of cliiraf mixing [35] on isovector- vector {V) and -axialvector [A) 
spectral functions. When extrapolated to the chiral restoration point (e=l/2), 
the degenerate V and A spectral functions also degenerate with the perturbative 
continuum (dashed line) down to s ~ 1 GeV^ , thus filling the "dip- region" in the 
vacuum V spectral function. 

onance. Such processes are of course included in the many-body treatment 
of the low-mass region discussed above, e.g., via pvr — )■ ai interactions. Here 
we use the mixing effect to predict that the dip region in the vector spectral 
function will be enhanced due to the admixture of the axialvector reso- 
nance (ai). In fact, if one takes the r decay data for the vacuum V and 
A spectral functions and carries the mixing all the way to the degeneracy 
point (e=l/2), one finds the intriguing result that the dip region is filled 
"precisely" to the level of the perturbative continuum, cf. Fig. [6l Stated 
differently, chiral mixing induces the reduction of the in-medium duality 
threshold to M^uai — 1 GeV. The mixing effect is not enough to smear the 
p resonance ~ this is where the (higher-order) hadronic many-body effects 
come in, as we will see below. Similar arguments can be made for nuclear 
matter at zero temperature, where the role of the thermal pions is taken 
over by the virtual pion cloud of the nucleon [2^ ; this is closely related 
to the pion cloud effects on the p meson through vrA^ scattering discussed 
above. 



2.3. Quark- Gluon Plasma and Lattice QCD 



In a quark-gluon plasma (QGP) the leading source of thermal dileptons 
is the purely electromagnetic annihilation of two quarks of equal flavor. 
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qq — >• e~^e~, corresponding to a structureless noninteracting EM spectral 
function, ImllEM oc M^, as given by the lower line in Eq. A consistent 
(gauge- invariant) calculation of loop corrections requires a careful treatment 
of infrared physics at the scale gT (g = ^/4^Ta^) which can contribute at the 
same order as the leading term (due to soft propagators of the type l/(t — 
m|,) with m£) ~ gT, canceling the extra coupling constants in the vertices). 
A systematic way to resolve these problems is provided by resummations 
within hard-thermal-loop (HTL) perturbation theory, which has been first 
applied to thermal dilepton production in Ref. [25J. The softening of the 
in-medium quark dispersion relation, developing plasmino branches as well 
as Landau damping, in the quark-antiquark loop produce an appreciable 
enhancement of the HTL dilepton rate (or EM spectral function) over the 
free one, not unlike the medium effects in the pion cloud of the p. The 
enhancement is present for a rather wide regime around qo ~ gT, and 
growing with decreasing go- 

In more recent years, the problem of soft dilepton emission has been 
addressed from first principles using lattice QCD (IQCD) [26], by computing 
the thermal expectation value of current-current correlation in Eq. ([3|) for 
Euclidean time, r, as 



In a heat bath in equilibrium, the imaginary-time coordinate, r, is not re- 
lated to dynamical but statistical properties of the system: the temperature 
comes in through a periodicity in r, effectively limiting it to the "Matsubara 
circle", < r < 1/T. The connection to the timelike (observable) spectral 
function can be established by means of a Kramers-Kronig dispersion rela- 
tion leading to 



where p{qQ,q;T) = —2 Imll^ is the polarization-summed spectral function 
(in the vector channel, charge conservation implies the 00-component to 
proportional to 5{qo)). Equation (fT2]l illustrates a notorious problem of ex- 
tracting spectral functions from Euclidean correlators computed in IQCD: 
one needs to perform an inverse integral transform on a finite number of 
meshpoints in r, which additionally carry statistical and systematic uncer- 
tainty and are limited to the Matsubara circle. Nevertheless, probabilistic 
methods (known as maximum entropy method, successfully used, e.g., in 
image reconstruction) can be applied which are usually coupled with addi- 
tional physical constraints such as positivity of the spectral function [27| . 




(11) 



oo 




(12) 
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Fig. 7. Left panel: EM spectral function, at vanishing S-momentum and summed 
over polarization, from thermal IQCD at T—lAbTc (solid black line) [2&, as well 
as for vacuum and in- medium hadronic calculations at T= 180 MeV (dashed and 
red solid line, respectively) [SB]. The spectral functions are rendered dimensionless 
upon division by energy, temperature and a charge/isospin degeneracy (Cem — 
Tliu d'^q ~ '^'-'^ IQCD and Cem — 1/2 for the isovector hadronic results). Right 
panel: euclidean correlators corresponding to the spectral functions in the left 
panel (squares: IQCD; solid lines: in-medium hadronic with the free continuum 
(lower line) and one with a threshold lowered by 0.3 GeV (upper line); dashed line: 
vacuum; the hadronic calculations additionally include the contribution from free 
or in-medium io spectral functions). 



However, it currently remains unclear as to how much structure can be 
resolved in the spectral functions (e.g., resonance peaks vs. threshold en- 
hancements). In Ref. [26] , this problem has been mitigated by making a 
physically motivated ansatz for the EM spectral function and fitting its 3 
free parameters to the euclidean correlator "data". The ansatz consists of 
a Breit-Wigner "transport" peak at low energies plus a pertubative contin- 
uum with as correction (plus Pauli-blocking factor) at high energies, as in 
the lower line of Eq. ([6]). The resulting spectral function in a gluon plasma 
of temperature 1.45 Tc (without thermal anti-/quarks) is shown as the solid 
black line in left panel of Fig. [7] It is divided by the energy variable to 
exhibit the "transport peak" for qq ^ which is dictated by the proper 
low-energy limit of the retarded correlator, PEuilo — )■ 0;T) oc crEuiT) qo, 
where ctem denotes the electrical conductivity (a transport coefficient of the 
strongly interacting medium). Also shown are the vacuum correlator and 
the results of the hadronic many-body calculations predicting a strongly 
broadened p meson at a temperature close to Tc. Despite the seemingly 
quite different environments (gluon plasma vs. interacting hadron gas) the 
two in-medium calculations are very different from the vacuum shape, but 
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not very different from each other. In fact, a cleaner comparison can be 
done by going back to the euchdean correlators, by simply performing the 
integral over the hadronic many-body spectral function in Eq. ()12p . The re- 
sulting correlators (normalized to the non-interacting one) for the vacuum 
and in-medium spectral function are compared to the direct IQCD com- 
putations in a gluon plasma in the right panel of Fig. [71 The in-medium 
spectral function leads to a correlator ratio quite similar to the one from 
IQCD, while the lack of a low-mass enhancement in the vacuum spectral 
function entails a ca. 50% underestimate of the maximum in n(r = 1/2T) 
(the euclidean correlators are symmetric about the midpoint of the Mat- 
subara circle, r = 1/2T). The rather good agreement of the in- medium 
correlators reiterates that the latter are not particularly sensitive to de- 
tailed structures in the spectral function, but that they do contain valuable 
information, such as a redistribution of strength to low energies or estimates 
for the transport peak. 



2.4- In-Medium Diletpon and Photon Rates 

It is now straightforward to convert the EM spectral functions dis- 
cussed in the previous two sections into thermal emission rates of dilep- 
tons, cf. Fig. [HI The left panel confirms that the strong broadening of the 
p spectral function, together with the chiral mixing in the dip region, make 
the hadronic rate approach the partonic-based calculations, in particular the 
HTL-improved result. In Ref. [28] this has been interpreted as an in- medium 
reduction of the quark-hadron duality in hadronic matter with increasing 
temperature and density: the hadronic rate successively approaches the per- 
turbative partonic rate, with Mjuai ~^ in the vicinity of Tc (or somewhat 
above). On the one hand, this implies an approach to chiral restoration, 
since the perturbative rate is chirally restored, i.e., degenerate with the ax- 
ialvector. On the other hand, it also implies an approach to deconfinement, 
as signaled by the transition to a more economic (perturbative) description 
in terms of partonic degrees of freedom. This picture is now further cor- 
roborated by the most recent IQCD results |26j : the right panel of Fig. [8] 
exhibits a rather close agreement of the nonperturbative IQCD rate with the 
HTL rate (except at very small energies go ^ gT where the latter receives 
additional corrections). 

In Ref. |29] , a "more realistic way" to interpret quark-hadron duality has 
been suggested, by advocating the 3-volume multiplied production rates, 
VdR/dM'^ , as the relevant quantities to compare. This is to account for the 
potentially rather different 3- volumes of hadronic and QGP phases if consid- 
ering a system at fixed entropy (which is a good approximation for a fireball 
in URHICs). However, the thermal production rates are thermodynamically 
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Fig. 8. Thermal dilepton rates; left panel: from a hadron gas with vacuum (solid 
green line) and in-medium [2 8) (solid red line) EM spectral function, compared to 
qq annihilation in leading-order (dashed line) and with hard-thermal loop correc- 
tions [25] (dashed-dotted line); all rates are integrated over the pair 3- momentum; 
right panel: from quenched lattice QCD at 1.45 Tc ^6] as extracted from the eu- 
clidean correlator shown in Fig. [7] (black and green lines), also compared to LO 
and HTL calculations, all at vanishing pair 3-momentum. Right figure taken from 
Ref. dS]. 



intensive quantities, well-defined in the infinite- volume limit (in which case 
one should not multiply with a volume factor). Furthermore, as discussed 
in the introduction, the finite-T QCD phase transition is most likely not a 
sharp one, especially concerning deconfinement. Recall that the increase in 
the degrees of freedom observed in the thermodynamic state variables such 
as energy density is well accounted for by the increasing population of states 
in the hadronic resonances gas [6j. It is precisely this increase (a precursor 
of the so-called the "Hagedorn" catastrophe) which leads to a likewise in- 
crease in the p width thus signaling its melting. In other words, the increase 
in degrees of freedom, which is ultimately responsible for a reduced volume 
at fixed total entropy in a finite system, is a hadronic effect, which, in turn, 
is at the origin of the p melting (as can be seen by glancing at the right 
panel of Fig. [TBI below). 
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Fig. 9. Thermal photon rates from hadronic matter |30) (sohd blue line: vector 
spectral function [28 at the photon point; dashed line: f-channel meson-exchange 
reactions) and from the QGP (solid red line: LO HTL calculation [21], dashed- 
dotted line: complete resummed LO result [5^). 

Thermal photon rates from the hadronic and QGP phase are compiled 
in Fig. [9l As elaborated following Eq. photon production processes 
from a medium of on-shell particles require a scattering process and are 
therefore nontrivial to leading order in the strong coupling (e.g., irp — )• tt'j 
or qg — t- q'j). However, in consistently resummed many-body calculations 
of the EM spectral function, the latter possesses a nonzero M — )• limit 
at finite 3- momentum, see, e.g., the p spectral function in the left panel of 
Fig. [5] (if additionally g — )• one obtains the electric conductivity) . This has 
been used in Ref. [3D] to extract the photon production rates encoded in the 
hadronic spectral function of Ref. |2H], corresponding to processes of type 
TT/O — )• ai — )• 7, Np — )• A^* — )• A^7 or pion t-channel exchange in vrA^ — )• 
included in T,pM, ^pB and SpTTTr, respectively (cf. the solid blue line in 
Fig. [9]). Additional t-channel meson-exchange reactions become important 
at photon momenta q > 1 GeV (cf. the dashed line in Fig. [9|), most notably oj 
exchange in np — )• ttj. These rates are compared to QGP-based calculations, 
specifically to the LO HTL rate derived in Ref. [3T] . 



for Nf=3 light flavors and with an added "1" in the logarithm to regularize 
its infrared behavior (surprisingly, the such obtained electric conductivity 




(13) 
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is not far from the one computed in thermal IQCD discussed in Sec. I2.3p . 
This expression nicely illustrates the leading 0{asCi) behavior of the photon 
rate. It turns out, however, that there are additional contributions at this 
order, requiring a full resummation of ladder diagrams. This has been 
achieved in Ref. ^32j , The numerical result of this work exceeds the "naive" 
LO rate by about a factor of 2-3, cf. Fig. [9j When comparing QGP and 
total hadronic rates (the sum of the two contributions plotted in Fig. [U]) , one 
again finds that both are very comparable at a temperature of r=200 MeV 
(chosen mostly for historic reasons). Thus photon rates also support the 
duality hypothesis, although the theoretical control over the calculations is 
less good as in the dilepton case. 

3. Dilepton and Photon Spectra in Heavy-Ion Collisions 

To date, dilepton spectra in heavy-ion collisions have been measured at 
the SPS (^=17.3,8.7 AGeV) [MIIMIEHIIMIEZIEH] at RHIC (\/s=200 AGeV) 
[391 SO] and at lower energies at BEVALAC and SIS {Ei^y,=l-2 AGeV) [HI 
142] . see, e.g., Refs. |43l l44j for recent experimental reviews. In this section 
we apply the rates discussed above to calculate thermal-emission spectra 
from URHICs, focusing our comparison on recent highlights from data at 
SPS and RHIC. Two additional ingredients are required to do so: a realistic 
space-time evolution of the expanding fireball formed in the reaction and 
an assessment of non-thermal sources contributing to the thermal spectra 
(in some cases it is possible to subtract these experimentally). 

Basic elements of a realistic fireball are summarized in Sec. 13. H while 
non-thermal sources are addressed on a case-by-case bases of each observ- 
able. The analysis of experimental spectra is organized according to their 
putative information content, namely in-medium effects on spectral prop- 
erties in connection with the emission duration (Sec. 13. 2p and temperature 
and expansion characteristics of the source (Sec. 13.3]) . 

3.1. Space-Time Evolution of Heavy-Ion Collisions 

Much has been learned about the space-time evolution of the medium 
formed in URHICs from the wealth of available hadron abundances, pt spec- 
tra and elliptic flow. Generally, hydrodynamic models do well in describing 
the latter two (especially at RHIC and LHC) up to momenta of pt=2-3 GeV. 
This encompasses more than 90% of the observed particles, and justifles 
the notion of a collectively expanding medium in local thermal equilibrium 
where the pressure gradients drive the acceleration. This mechanism works 
well enough to have predicted azimuthal asymmetries in the pt spectra for 
non-central heavy- ion collisions: the elliptic asymmetry given by the trans- 
verse overlap zone of the colliding nuclei leads to larger pressure gradients 
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Fig. 10. Schematic representation of isentropic trajectories that the medium fol- 
lows in URHICs at SPS, RHIC and LHC in the QCD phase diagram (left panel) 
and when projected on the temperature axis as a function of time using hadron 
spectra to estimate the collective expansion properties (right panel). Each trajec- 
tory is calculated for a fixed specific entropy, S/Nb — s/qb, according to entropy 
and baryon-number conservation. In addition, the hadronic part of the trajecto- 
ries implements effective chemical potentials for all strongly stable hadrons [33] 
to maintain the observed hadron ratios as determined by thermal-model fits at 
chemical freezeout [I] characterized by the "data points" . 

along the short axis of the ellipsoid, which in turn provide larger accelera- 
tion resulting in a positive "elliptic flow" coefficient, V2{pt), in the angular 
distribution of the hadron pt spectra (this is arguably the strongest evi- 
dence for the celebrated "near-perfect liquid" discovered at RHIC). Taking 
all evidence together, the following picture emerges: Shortly after the col- 
lision of the two nuclei, at a thermalization time of to=0. 5-1 fm/c, most 
of the entropy has been produced and the system has thermalized into a 
QGP, with initial temperature of ca. ro=350-400 MeV (200-250 MeV at SPS, 
~500-600 MeV at LHC). The QGP expands for about 3 fm/c before entering 
the transition regime to hadronic matter at around Tc — 170 MeV, where 
it spends almost 5 fm/c. At this point the hadronic abundances become 
frozen, leading to the notion of "chemical freezeout" where the inelastic re- 
actions are believed to turn off. In the subsequent hadronic phase, elastic 
scattering driven by hadronic resonance is still operative leading to a further 
cooling of the system for about 5-10 fm/c until "thermal freezeout" occurs 
at Tfo ~ 100 MeV. After this, hadrons stream freely to the detector (modulo 
long-lived electromagnetic and weak decays). Figure [TO] summarizes typical 
trajectories in the phase diagram, as well as the time dependence of tem- 
perature, for central heavy-ion collisions at SPS, RHIC and LHC. Besides 
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initial conditions and the equation of state, hydrodynamic models require 
the implementation of "freezeout" . At thermal freezeout, the medium is con- 
verted into hadron spectra using the time-honored Cooper-Fry prescription. 
"Chemical freezeout" implies that in the subsequent evolution the abun- 
dances of observed (stable) hadrons should be conserved. This is achieved 
by introducing effective chemical potentials into the hadronic equation of 
state (usually modeled via a hadron resonance gas), to individually con- 
serve the number of pions, kaons, nucleons and antinucleons (not only their 
difference, which is an exact conservation law related to the baryon chemi- 
cal potential, /x^), etc. [l5j. For dilepton (photon) production this implies 
that, e.g., the process vrvr — )• p — )• e+e~ (vrp — )• -k^) is augmented by a pion 
fugacity squared, {z^.) where z-j^ = exp(^^/T). 

Before going into quantitative theory comparisons to data, it is instruc- 
tive to make a more general assessment of which stages of the fireball evo- 
lution will contribute at which dilepton energies / masses in the spectra. To- 
ward this goal we write the the time-differential invariant-mass spectrum, 
integrated over 3-momentum and 3- volume at each time snapshot, as |10] 

C^ATee _ M r . 3 diVee 

dMdT ~ qoJ '^ d^xd^q 

^ ' M J Qo 

^ const VMr) ^"^"^^^f^f'^^ e"^^/^ (MTf/^ , (14) 

where VpB denotes the volume of the expanding fireball. For the second 
(approximate) equality it has been assumed that the in-medium EM spectral 
function depends only weakly on 3-momentum (which is exact in vacuum 
due to Lorentz-invariance, and approximately satisfied in medium for not 
too large momenta). For the third (approximate) equality we have invoked 
the non-relativistic approximation, T/M <C 1. In a last step, we relate the 
fireball-evolution time to the temperature, T{t), to obtain the temperature- 
differential emission as 

"^^^'^ ocImnEM(M;r) e-*^/^ r— . (15) 



dMdT 



The power of temperature, m, depends on the details of the fireball ex- 
pansion (e.g., 1-D Bjorken or full 3-D as appropriate for the early and 
late phases, respectively), and the equation of state (i.e., the power law of 
entropy density, s oc T", with n=3 for an ideal massless gas, but larger 
for a hadron-resonance gas), see Ref. [10] for more details. In Fig. [TT] we 
plot the coefficient /(T;M) = exp(— M/T) T~"^ as a function of temper- 
ature for different masses. At each mass, the maximal emission yields are 
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Temperature Profile of Dilepton Yields 
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Fig. 11. Emission yields of thermal dileptons as a function of temperature for 
different values of their invariant mass, assuming a temperature independent EM 
spectral function. 

centered around a temperature T^ax — T/m ~ T/S which reflects the com- 
petition between the thermal exponential, favoring high temperature, and 
the increasing 3-volume of the fireball as it expands and cools. Note that 
for low-mass dilepton emission, M < 1 GeV, the profile function, f{T;M), 
peaks in the hadronic phase of the evolution. On the one hand, medium 
effects in the EM spectral function, which lead to a marked enhancement 
for M < rup with increasing T and qb (recall Fig. [5]), will shift this peak to 
somewhat larger T (note, however, that the low-mass enhancement tends 
to saturate when approaching Tc, see Fig. [8] left). On the other hand, the 
latent heat burned off in the hadron-to-parton transition, i.e., the rapid rise 
of s/T^ in the vicinity of the phase conversion, will bias the emission to 
the hadronic phase. Both effects lead to an additional focusing of low-mass 
dilepton emission to hot and dense hadronic matter. 

Similar arguments apply to thermal photon (or dilepton) transverse- 
momentum spectra, with an important amendment due to the transverse 
fiow of the fireball causing a spectral blue shift (well-known from hadron 
spectra). Since the flow develops over time, it creates a bias toward the 
later emission stages. 

3.2. Spectrometer and Chronometer 

The state-of-the art in dilepton measurements in URHICs has been set 
by the NA60 collaboration with their dimuon data from In-In(y^=17.3 AGeV) 
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Fig. 12. Left panels: excess dimuon invariant-mass spectra as measured in In- 
In(Y^=17.3 AGeV) collisions at the SPS [23 EH], compared to theory calculations 
with in-medium hadronic spectral ftmctions and QGP emission from an expanding 
fireball [13]. The upper panel resolves the individual contributions in the calcu- 
lation, which, however are dominated by the in-medium p spectral fiinction at 
low-mass and a combination of hadronic continuum ("47r") and QGP radiation at 
intermediate mass. The in-medium p plus 47r make up the total in-medium hadron- 
gas (HG) rate shown by the solid red line in the right panels for two temperatures 
(including baryon-density effects) [TS] . The lower left panel illustrates the effect of 
varying the fireball lifetime in the calculation. 



collisions at the SPS |35l [37]. Excellent statistics combined with superb 
mass resolution allowed for a subtraction of background sources thus isolat- 
ing the "excess radiation" from reinteractions in the fireball. In addition, 
a full acceptance correction could be carried out so that the invariant-mass 
spectra, displayed in the left panels of Fig. [121 are, for the first time, truly 
invariant. This, in particular, implies that blue-shift effects of the expand- 
ing fireball do not distort the spectra. The result is quite stunning: the data 
show a near-perfect thermal spectrum, "only" modulated by a broad bump 
around the free p/uj masses. Even without any fireball convolution, a com- 
parison to the theoretical input rates calculated in Ref. [12] (shown in the 
left panel of Fig. [SI supplemented with dimuon threshold kinematics) ex- 
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Fig. 13. Left panel: comparison of the acceptance-corrected NA60 mass spectra in 
In-ln(yi=17.3 AGeV) collisions [371 [35] to a calculation where baryon-induced ef- 
fects ion the EM spectral function have been switched off. Right panel: in-medium 
p widths (estimated from the full-width-half-maximum of its spectral function) 
along the SPS trajectory used in the spectra calculations projected on the temper- 
ature axis, for the full calculation including effective chemical potentials (circles) 
and in chemical equilibrium with (squares) and without (triangles) baryon-induced 
effects; figure taken from Ref. (4]. 

hibits a remarkable agreement for an average temperature of 150-160 MeV. 
This temperature is tantalizingly close to the most recent IQCD predictions 
shown in the left panel of Fig. O provoking the conclusion that the chiral 
transition has been observed at SPS! Of course, these rather general argu- 
ments need to be backed up by quantitative calculations convoluting the 
T- and ^^-dependent emission rates over a realistic fireball expansion (as 
discussed in Sec. 13. ip . This has been done in Ref. 06] and the results are 
in very good agreement with the data as seen in the upper left panel of 
Fig. [12j Subleading contributions include the decays of uj and (f) mesons, 
a non-thermal component of p mesons, and Drell-Yan annihilation. Espe- 
cially the latter two play a significant role when high-gj cuts are imposed 
on the spectra (this is fully in line with the standard notion from hadron 
Pt spectra that the thermal component dominates up to pt —2 GeV yielding 
to "hard production" above). With the level of precision set by the NA60 
mass spectra several further diagnostics of the produced medium become 
available, as discussed in the following. 

The first obvious question is how sensitive the mass spectra are to the p- 
meson spectral shape. This has been tested qualitatively by switching off the 
medium effects induced by baryons, i.e., using the T-dependent p spectral 
function in a hot meson gas as displayed in the right panel of Fig. [SJ The re- 
sulting p^p~ spectra are in stark disagreement with the NA60 data, cf. left 
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panel of Fig. [T3l the low-mass enhancement is underestimated by about a 
factor of ~4, and the yield in the p-peak region is overestimated by about a 
factor of ~2, both several standard deviations away from experiment. The 
evolution of the p width (taken as full-width-half-maximum in the spectral 
function) along the SPS trajectory used in the calculations is displayed in 
the right panel of Fig. [13l The upper curve (circles), corresponding to the 
full calculation, confirms that the "average width" of V^'^'^ ~ 350-400 MeV 
reflects the conditions of the medium close to the chiral transition. It also 
exhibits a rather strong rise for temperatures beyond that. This rise be- 
comes even more pronounced in chemically equilibrated matter where the 
effective meson-chemical potentials are switched off (squares; at fixed s/qb 
this also implies a slight reduction va. ^b)- It is not unlike the increase of 
the thermodynamic state variables in the hadron-resonance gas around a 
similar temperature. No such behavior emerges at comparable tempera- 
tures if baryon-induced medium effects are switched off. Since the average 
width implies the occurrence of larger widths in the course of the medium 
expansion, a situation with F™^'^ ~ nip must be realized, i.e., resonance 
melting. 

It is important to note that baryon-induced medium effects remain a key 
player at RHIC [47J . This is so because the p equally interacts with baryons 
[B) and anti-baryons {B) due to CP invariance of strong interactions (the 
p is a CP eigenstate). Therefore, the sum (not the difference) oi qb and 

is the relevant quantity for medium effects. At RHIC, the rapidity 
density of baryons plus antibaryons is slightly larger than at SPS, while the 
pion rapidity density is roughly a factor of 2 larger, implying an accordingly 
increased fireball volume at given temperature. Thus, the "effective" density 
Peflf = Qb + Qb-> is very different at SPS and RHIC, and probably also 
at LHC. From a theoretical point of view, the configuration of small net- 
baryon chemical potential at collider energies is appealing because it comes 
closest to pb=0 where IQCD computations are most powerful. 

While the spectral shape of the NA60 data could be predicted by the- 
ory, the total yields could not, at least not within the experimental accu- 
racy. However, one can turn the argument around and utilize the total 
yield as a precision clock for the emission duration. In previous analyses 
of dilepton spectra in Pb-Au collisions by CERES [33], the fireball lifetime 
estimated from hydrodynamic models, rpB — 10 — 15fm/c, resulted in a 
reasonable description of the e^e~ excess radiation p8j . The precision at- 
tained with the NA60 data is significantly increased: even a variation of 
=blfm/c in the calculation discussed above leads to noticeable deviations 
from the data, cf. lower left panel of Fig. [121 The resulting "measurement" 
of TrrriFB = 6.5ib Ifm/c for the fireball lifetime using dilepton yields is prob- 
ably the most accurate one thus far in URHICs (it is not to be confused 
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with the pion emission duration referred to in Hanbury-Brown Twiss (HBT) 
inter ferometry measurements). 

Let us briefly comment on the current situation of low-mass dilepton 
measurements at RHIC. The PHENIX collaboration found a very large en- 
hancement for masses below the free p mass in central Au-Au(-y/s=200 AGeV) 
collisions [39] which is largely concentrated at low transverse momenta, 
qt < 0.5 GeV. It cannot be explained by the theoretical ingredients which 
allow for an understanding of the SPS dileptons. The STAR collaboration, 
on the other hand, reports a smaller enhancement [iQ], which, however, is 
still sizable and roughly consistent with theoretical expectations. 

3.3. Thermometer and Barometer 

Historically, the best experimental tool to serve as a "thermometer" of 
the early phases in URHICs was believed to be thermal radiation of photons. 
With a structureless emission rate, the slope of the spectrum is directly re- 
lated to the temperature of the system, modulo a blue-shift from transverse 
expansion which is expected to be small in the early stages. At SPS en- 
ergies, the presence of a significant excess signal of photons, i.e., beyond 
those from initial production and final-state decays of long-lived hadrons 
(mostly TT^jT] — )• 77), remains inconclusive [48j. At RHIC, however, the 
PHENIX collaboration has observed a strong excess signal in semi-/central 
Au-Au collisions [49j for transverse momenta above qt = 1 GeV up to at 
least 3-4 GeV, beyond which the primordial production extrapolated from 
binary A^-A^ collisions dominates, see left panel of Fig. [T3] ("excess" plus 
"primordial" are usually referred to as "direct" photons, i.e. long-lived de- 
cays subtracted). This range nicely coincides with theoretical expectations 
of where the QGP radiation is most visible [M EQl EQl ED IMl |53] . The 
inverse slope of the excess radiation was found to be consistent with an 
exponential with temperature T = 221 it 19^*^* ± I9^y^ MeV. This seems to 
point at a source well inside the QGP, albeit not to the very early phases 
where temperatures above 300 MeV are expected, recall Fig. [TOl 

Maybe even more intriguing are the recent PHENIX measurements of 
the direct-photon elliptic fiow [5ll , cf. right panel of Fig. [TH In the regime of 
the excess radiation, i.e., for qt < 3 GeV, the elliptic-flow coefficient, V2{qt), 
is comparable to those of pions, which are only emitted at the end of the fire- 
ball evolution, i.e., with maximal V2. On the contrary, the excess photons are 
supposed to be emitted throughout the evolution (including early, where the 
V2 is still small) , and the primordial ones should carry zero f 2 . More quanti- 
tatively, predictions of thermal-photon f 2 using hydrodynamic models with 
QGP-dominated excess radiation lead to maximal values of V2 — 2% for di- 
rect photons in 0-20% Au-Au(Vs=200 AGeV) collisions [Ml Ell E3 153] , weU 
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Fig. 14. Transverse-momentum spectra (left panel) and elliptic flow (right panel) 
of direct photons in 0-20% central Au-Au(-ys=200 AGeV) collisions at RHIC, com- 
pared to PHENIX data [IHl IM] • The curves are calculations [5^ with a realistic 
fireball evolution employing thermal QGP and hadronic rates which are "dual" 
around Tc, corresponding to Fig. El 



below the PHENIX data. In these calculations the hadronic emission rates 
were significantly smaller than the QGP rate in the vicinity of Tc, which 
differs from the "duality" hypothesis suggested by Fig. El In Ref. [55J we 
therefore revisited our earlier calculations of direct photons ^30j by refining 
the thermal fireball model as to quantitatively describe the pt spectra and el- 
liptic flow of bulk hadrons, freezing out at Tfo — 100 MeV, and multistrange 
hadrons, kinetically decoupling close to chemical freezeout, Teh — 180 MeV. 
This, in particular, led to harder thermal photon spectra from the hadronic 
phase, thus exceeding the QGP contribution over the for thermal emis- 
sion relevant qt range up to 3-4 GeV. Consequently, the thermal-photon V2 
increases substantially, now being dominated by hadronic emission where 
most of the final V2 has already been built up. The resulting direct-photon 
^2 (<?*)) which includes the primordial component (extrapolated from pp col- 
lisions) with vanishing V2, reaches up to 6-7%, a factor of ~3 larger than in 
previous calculations and thus much closer to experiment, cf. right panel in 
Fig.m 

With dominantly hadronic emission, one should revisit the evaluation of 
the photon slope parameter. For this purpose, we plot in the left panel of 
Fig. [151 the blue-shifted effective temperature as a function of time in the 
underlying fireball evolution. 



where T is the "true" temperature in the thermal rest frame and (/?) the 
average transverse flow velocity (typically 2/3 of the surface velocity). The 
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Fig. 15. Effective slope parameters of transverse-momentum spectra of EM ra- 
diation in URHICs. Left panel: time dependence of T^s of thermal photons, 
Eq. in an expanding fireball for 0-20% central Au-Au(ys=200 AGeV) colli- 

sions at RHIC (solid blue line), compared the local temperature (dashed line) and 
PHENIX data [S3] (horizontal solid line with dashed-line error band); the dash- 
dotted line indicated the flow velocity of the fireball surface. Right panel: Toff for 
excess dimuons as a function of their invariant mass as measured by NA60 in In- 
ln(Y^=17.3 AGeV) collisions at SPS [37], compared to theory calculations [3S] with 
different transverse fireball acceleration (lower 2 curves: at =0.085 fm/c^, upper 2 
curves: at =0.1fm/c^) and with different phase-transition temperatures (EoS-B: 
Tc^ieOMeV, EoS-C: Tc^igOMeV). 

comparison to the PHENIX measurement reveals that the largest overlap is 
indeed inside the hadronic phase, for a broad temperature interval around 
~130MeV, "Doppler-shifted" by about 100 MeV. This is consistent with 
the large elliptic flow. 

Similar analysis can be carried out for dileptons, which has been done 
with the NA60 dimuon gj-spectra at the SPS [37j. The extracted slope pa- 
rameters are displayed as a function of the invariant-mass window imposed 
on the qt spectra in the right panel of Fig. [15] One finds a continuous rise of 
Tefj with mass, roughly following the schematic (nonrelativistic) analogue 
of Eq. ()16p for massive particles (here: virtual photons), 

reff~r + M(/3)2 . (17) 

This corroborates the hadronic nature of low-mass dileptons below the p 
mass. However, above the p mass, the slopes decrease, which has been 
argued to signal a transition to a new source with small flow and thus of 
partonic nature. However, the latter part of the interpretation is not unique 
and one should keep in mind the following features: (a) the p-mass region, 
M ~ 0.7-0.9 GeV, is predominantly populated by contributions from the 
late flreball stages, since the p spectral shape is recovering its vacuum form 
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which prominently figures around M ~ nip] (b) the interplay of thermal fac- 
tor and 3- volume implies that going up in mass the emission is "biased" to- 
ward the earlier phases (recall Fig. llOp . which increases T but more strongly 
decreases the Doppler shift in Eq. (jl7p . M{/3)'^; (c) for temperatures around 
Tc the duality of the emission rates renders a "partonic" and "hadronic" 
assignment ambiguous, depending on the not-so-well-defined definition of 
the transition temperature. The latter is demonstrated by the two pair of 
lines in Fig. [TC} where each pair is calculated assuming two "extreme" val- 
ues of rc=160 and 190 MeV. In the former case, the total thermal yield for 
masses M=l-1.5 GeV is indeed dominated by QGP emission, while in the 
latter case hadronic emission prevails (with chiral mixing, cf. Fig. [6]). The 
corresponding two curves in a pair essentially overlap. A more significant 
effect is due to transverse flow. The original predictions of fireball [57] 
and hydrodynamic [ 58] evolutions agree with each other, but underestimate 
the data \59\ (cf. the lower pair of curves in Fig. [15] right). When imple- 
menting a stronger transverse expansion [36], comparable to that used at 
RHIC {at=0.1 fm/c^), a reasonable description of the NA60 slopes can be 
obtained (upper two curves in Fig. 1151 right). Using the schematic slope for- 
mula with T=175MeV and an average QGP+mixed phase evolution time 
of 3fm/c for semicentral In-In (resulting in (/3) ~ |a(T = 0.2c) gives an 
inverse slope of Tefj ~ 215 MeV for M ~ 1 GeV, as seen in the figure. The 
main point is that the NA60 slopes for M > 1 GeV give an independent 
confirmation of thermal emission from a QCD medium with dual rate and 
temperature of around Tc. The sensitivity to the unexpectedly large col- 
lective fiow shows that dileptons can also serve as an accurate "barometer" 
of the fireball. A similar analysis at RHIC should lead to much increased 
slope parameters, since the QGP has built up much more fiow by the time 
it has cooled down to Tc. For example, from the left panel of Fig. [15] 
one finds f3s — 0.6fm/c^ at the end of the mixed phase, translating into 
Tefj ~ 350 MeV, see also Refs. [60\ 161]. Finally we remark that invariant- 
mass (rather than transverse-momentum) dilepton spectra can, in principle, 
provide a cleaner measurement of the temperature since it is Doppler-free. 
Of course, the richer (dynamic) structure in the low-mass regime complicates 
this task, but with a reliable knowledge of the in-medium spectral shape it 
is still possible; for the structureless regime above M=l GeV the situation 
is even better. Analyzing the invariant-mass slopes of the NA60 spectra 
one finds temperatures ranging from 150-180 MeV for M ~ 0.2 — 1.5 GeV, 
approaching 200 MeV for M > 2 GeV. The latter regime requires a good 
knowledge of the primordial contribution, i.e., Drell-Yan annihilation and 
its nuclear corrections 1621. 
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4. Conclusions 

Electromagnetic emission spectra provide a rich observable to investi- 
gate the strongly interacting medium produced in high-energy collisions of 
heavy nuclei. The starting point for exploiting this observable is a thorough 
understanding of the vector spectral function in equilibrium matter. Chiral 
effective models at low and intermediate temperature and densities, lattice- 
QCD in the vicinity of the phase transition(s) and perturbative QCD in the 
high-tcmpcraturc limit are starting to fill the theoretical landscape. The 
evidence is mounting for a melting of the hadronic vector-meson resonances 
and a transition to a partonic-like rate over the entire mass range. This, 
in particular, includes a large low-mass enhancement which is corroborated 
by recent lattice QCD results for cuclidcan correlators and has been con- 
sistently showing in dilepton measurements from BEVALAC/SIS via SPS 
to RHIC. The state-of-the art in dilepton experimentation has been set by 
the NA60 dimuon results, which accurately confirm the resonance melting 
and show that it occurs right around the temperatures where IQCD predicts 
the chiral phase transition. This sets the stage for exciting developments 
to come, at both collider energies (RHIC and LHC, where the connections 
to IQCD are even more direct) and at lower energies as planned at the 
Compressed Baryonic Matter experiment (CBM, where the baryon-induced 
medium effects are expected to become max;imal). Several theoretical de- 
velopments will have to accompany these efforts, including hadronic calcu- 
lations of the axialvector to tighten the grip on chiral restoration using sum 
rules, lattice calculations of euclidean correlators with dynamical quarks, 
and the implementation of equilibrium rates into state-of-the-art hydrody- 
namic models (possibly supplemented with transport simulations) which 
quantitatively describe bulk-hadron observables. One can thus hope that 
the systematic use of EM probes will continue to deliver unique and precise 
information on the phase structure of a strongly coupled non-abelian gauge 
theory. 
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